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Abstract. We use the full multiple scattering expansion of the retarded self- 
energy to obtain the gain and loss rates present in the Kadanoff-Baym relativis- 
tic transport equation. The rates we obtain include processes with any number 
of particles. As a first approximation, we only consider central cuts in the 
self-energies, but otherwise our results are general. We specialize to the case of 
scalar field theory to compare with lowest order results. The main application 
of this work is relativistic transport theory of very dense systems, such as the 
quark-gluon plasma or the early universe, where multi-particle interactions are 
important. 
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1. Introduction 



Transport theory is an important tool in the study of the quark-gluon plasma (e.g. 
[IB El)- I n the case of the quark-gluon plasma (QGP), it is used to model some 
aspects of its evolution. In particular, in QGP thermalization studies, many authors 
(e.g. [0] |SJ H>]) stress the importance of including multi-particle interactions such as 
gg — > ggg in transport theory models. This is to be expected, since at the beginning 
of the QGP expansion, its energy density is > 1 GeV/fm 3 (approximately 6 times 
the density of cold nuclear matter). The same conclusion is drawn in transport co- 
efficients computations E]- Thus including multi-particle interactions in kinetic 
equations seems to be an important ingredient in the study of very dense systems. 
The goal of this paper is to show how to obtain complete relativistic gain/loss 
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rates containing any number of particles. These rates are the ones present in the 
collision term of a relativistic Kadanoff-Baym transport equation. Our derivation 
is based on the exact multiple scattering expansion of the self-energy [ |S] . 

2. Relativistic Transport Equations 

The results presented in this section rely on the work of [ 7 . In this paper, the 
authors derive a relativistic transport equations from quantum field theory for a 
system of real scalar fields (cubic and quartic interactions) : 



Bphr + V'V )f(X,p) + VV(X,p)V p f(X,p) = 6(p )C(X,p) (1) 



along with the mass-shell condition p 2 — m 2 + V(X, p) = 0, where E p is the solution 
of the mass-shell equation, v = dE p /dp is the velocity, f(X, p) is the distribu- 
tion function, V(X, p) is the Vlasov term (corresponding to an effective mass) and 
C(X,p) is the collision term (given below). To obtain Eq. ^ an d the mass-shell 
condtion, one must assume weak inhomogeneity with respect to the inverse charac- 
teristic momentum (allows for the gradient expansion) and the Compton wavelength 
(justifies the quasi-particle approximation) and also weak coupling. The collision 
term in Eq. \T]is given by: 



where U < ' > (X,p) are Wightman functions. These Wightman functions are "com- 
pletely cut" self-energies and are interpreted as gain/loss rates. These microscopic 
collision rates must be computed from quantum field theory. In [[7], the authors 
use perturbation theory and compute the rates explicitly. In this paper, we instead 
use the full Multiple Scattering Expansion of the Self-Energy to derive the rates 
and express them in terms of physical processes. 

3. The Multiple Scattering Expansion of the Self-Energy 

The idea behind the expansion is simple. For simplicity, let's consider only scalar 
theories. The starting point is the common structure of the four bosonic propagators 
corresponding to the contour Green's function in the real time formalism. Each 
propagator can be decomposed into a zero and finite temperature part, i.e. G %3 (p) = 
Gq (p) + T(p), the latter being common to the four propagators. Using this property, 
it is possible to reorganize any self-energy diagram as a power series in the number 
of thermal phase space factor T(p) — n(\p°\)2ir5(p 2 — to 2 ), where the coefficients of 
the expansion are zero temperature scattering diagrams 

To compute the rates n <:> (£;) appearing in Eq. [3 we need the multiple scat- 
tering expansion of the imaginary part of the retarded self-energy. Following the 




C(X,p) 



l -[u<(x,p)(f(x,p) + i)-n>(x,p)f(x,p)} (2) 
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Fig. 1. One-loop self-enegy example of multiple scattering expansion. The white 
circles correspond to thermal lines. The first and second terms are first order in 
the thermal expansion, while the third term is second order (the zero temperature 
term is not shown) 

procedure in [ [S] and concentrating on the loss rate, the thermal contribution to 
II > (fc) can be written as (a similar expression exists for n < (fc) [E]): 

00/ n+i „ \ 

n>(A) - c— E II / ^0 ( k dli}\tf%\k,{k}} C onn (3) 

where k is the external momentum, ij's are momenta corresponding to thermal 
lines, S 1 ^;} are symmetry factors, dTf — ^ 2 ^2E- n (-^)i Ei = Vh + m ii Ik 1S the 
part of the T matrix involving k (when it is in the initial state) and X)^°=o ^ s ^ ne 
summation over the expansion parameter dTj. The subscript "conn" means fully 
connected diagrams only. The operators inside the bra-kets are related to the usual 
T matrix containing all interactions. 



4. Gain/Loss Rates 

Let us consider the loss rate. The idea is to start from Eq. and insert a complete 
set of many-particle states between the two 7& operators. This gives: 



00 00 



n> w = EEMn 



nl 

n=0 s=0 




(2tt) 3 2Ei \ s! H J (2tt) 3 2E' 4 

/ \ j — 1 

x (k, {U}\^\{pj})({pj}\T k \k, {y)con„ (4) 

where the l/n\ comes from the symmetry factor of the cut diagrams It is 

understood that a product without a factor is one. This operation has the effect 
of "cutting" the diagrams on the RHS of Fig. ^ Let us concentrate on "central 
cuts" only, i.e. diagrams where only the lines that are already cut can be replaced 
by r(p); we treat the case of non-central cuts in another paper For the central 
cuts case, rif number of U's pass through Tk (they are spectators) and replace n/ 
of the Pj's. Inserting the factor 1 = ^n-n f -m^{ n > n f) to reorganize the sum 
and after a bit of algebra, we get: 



x X 2 

n ^cntral( fc ) = EE ^~T S ! 

s=0 ni =0 *' ' \i=l 

x(fc,{ii}iMi{Pi})<{Pi}i^fci*.w> 
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where the 7& matrix elements become A4k after all spectators are passed through 
(the S function associated with Mk is not written explicitly). Equation [S] is ex- 
pressed in terms of physical scattering processes and contains any number of par- 
ticles (the indices s and rij specify the number of final and initial particles, respec- 
tively). A similar expression can be obtained for the gain term [[§]. We can see 
that it reproduces the two- loop and three-loop results of [0 for n, ; = 1, s = 2 and 
rt,: = 1, s = 3. 

5. Conclusion 

We have shown that the multiple scattering expansion of the retarded self-energy 
can be used to obtain the gain/loss rates of the relativistic Kadanoff-Baym equation. 
These rates have a direct interpretation in terms of physical processes. Our approach 
has several advantages: the calculation is completely general, it contains any number 
of particles, it is valid to all orders and it is in principle generalizable to gauge 
theories. However, there is an important caveat: we assume that pinch singularities 
vanish, even if the system is out-of-equilibrium. This assumption is not true in 
general and we are currently working on improving this aspect of the problem. 
Details of the derivation and the generalization to non-central cuts will be discussed 
in another paper [Ej. 
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